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Abstract 

The regularized vacuum energy (or energy density) of a quantum field subjected to static external 
conditions is shown to satisfy a certain partial differential equation with respect to two variables, 
the mass and the "time" (ultraviolet cutoff parameter) . The equation is solved to provide integral 
expressions for the regularized energy (more precisely, the cylinder kernel) at positive mass in 
terms of that for zero mass. Alternatively, for fixed positive mass all coefficients in the short- 
time asymptotics of the regularized energy can be obtained recursively from the first nontrivial 
coefficient, which is the renormalized vacuum energy. 
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BACKGROUND 



The experimentally documented physical relevance of vacuum energy involves the elec- 
tromagnetic field (the Casimir effect). Other theoretical work also usually considers massless 
fields (scalar fields for calculational and conceptual simplicity; neutrinos or gluons for more 
exotic or speculative applications). Nevertheless, the vacuum energy of a quantum field 
with mass is of theoretical importance. Adding the mass term is the simplest generaliza- 
tion of the basic model of a massless scalar field in an empty fiat cavity, a step toward the 
more complicated scenarios with space-dependent external potentials or curved background 
geometries. Although massive Casimir-type effects in the simplest geometries tend to be 
exponentially darnped relative to their massless counterparts, they can be less negligible in 
other geometries Sees. 7.4-7.5]. 

The massive field in one-dimensional space was studied by Hays 0. (See also 0,30.) 
Perhaps the most notable feature of the results is the presence of a logarithmic divergence, 
absent from the massless case. Although physically harmless in the context of the one- 
dimensional bag model 0,0, a logarithmic divergence is regarded as particularly problematic 
by some theorists because it is not automatically eliminated by dimensional or zeta-function 
regularization. 

As I have stressed elsewhere 0, 0) vacuum energy is one of a series of moments of the 
spectral distribution of the differential operator appearing in the field equation; these quan- 
tities arise as the coefficients in the short-time asymptotics of a certain Green function, the 
"cylinder kernel". The diagonal value of the latter is, up to a differentiation and a numerical 
factor, the regularized vacuum energy defined by an ultraviolet cutoff (with the time vari- 
able as the cutoff parameter) 0, . Thus the vacuum energy (regardless of experimental 
relevance) is of mathematical interest as a tool of spectral analysis. The cylinder-kernel 
coefficients incorporate nonlocal geometrical information that is not extractible from the 
much-studied short-time asymptotic expansion of the heat kernel. (The cylinder kernel is 
minus twice the t-derivative of the Green function employed in 0,0.) 

NOTATION 

Consider a field equation of the type = —Hcj), where 

H = Hq + ^, /X = m^, (1) 

and Ho is a self-adjoint second-order differential operator in the spatial variables, such as 
— V^. For simplicity (although the treatment of local energy density is actually more general) 
assume that the spatial domain is compact and H has a discrete, positive spectrum {cu^} 
with orthonormal eigenf unctions {0„(x)}. 

For auxiliary mathematical purposes one studies the heat kernel 



K{t,x,y) = (x|e"*^|2/) = J] e-*'^"0.(x)0„(2/)* (2) 

n=l 

and the cylinder (Poisson) kernel 



T{t,x,y) = {x\e-'^\y) = ^ e-*-"0„(x)0„(y)*. (3) 

n=l 
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Each of these can be "traced" over space; for example, 

TrT = / (x|e-*^|x) dx = J2 e"*"""- (4) 

n=l 



oo 



Formally, the vacuum energy of the quantized field configuration is 

1 °° Id 
E= -V = --lim — TrT, (5) 

and one possible definition (see Hi) of the vacuum energy density is 

1 °° Id 

n=l 

In reality, the limits in (jSJ and © do not exist, but TrT and T(t, x, x) possess asymptotic 
expansions as t | of the form 0,13, HO, lOl 

oo oo 

T~^e,t-^+^+ fst-'^-'lnt, (7) 

s=0 s=a!+l 
s—d odd 

where the coefficients of the divergent terms are simple, local objects that can be absorbed 
by renormalization. (Here d is the spatial dimension.) Therefore, one regards Q and 
after operation by — , as the regularized energy and energy density, and one regards — | 
times the coefficient of the term of order t in Q and (j^ as the renormalized energy and 
energy density: 

E or Too = -|ed+i. (8) 
Similarly, if K stands for either Tr K or K{t, x, x), it has an expansion of the form 

oo 
s=0 



THE MAIN EQUATION 

The coefficients in ((7j) and Q are functions of /i. Let us write T(/i, t) and K{^, t) for the 
quantities being expanded and write es(/i), etc., for the coefficients. In the case of the heat 
kernel, it is elementary that 

if(/i,t) = ir(0,t)e-^*, (10) 

and from (jiup it is routine to find formulas for the &s(/i) in terms of 6s'(0) (s' < s). For the 
cylinder kernel it is clear that no elementary factorization like (|iup occurs, and hence the 
mass dependence is much more interesting and nontrivial. 

On the other hand, (jlUj) is equivalent to the differential equation 

dK , , 
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The goal of the present paper is to find, as nearly as possible, an analogue of (fTTjl for the 



quantities T(/i,t) related to vacuum energy. Since uJnifJ') = Y^n(O)^ + /i, it is easy to show 
from dH) or © that 



dfidtKtJ 2 ' 

which is the central equation of this paper. The variables t and n naturally range from to 
+00. 

If its right side were zero, (jl2|) would be mathematically equivalent to the massless wave 
equation in two-dimensional space-time written in null (light cone) coordinates; as is well 
known, its general solution would then be T{fi,t)/t = A{t) + B{fi), where A and B are 
arbitrary functions. The full equation (fT^ is of the same hyperbolic type, and one can 
again expect the general solution to involve two arbitrary one-variable functions. One of 
these should be the "initial value" T{0,t), in analogy with ()10j). The remaining boundary 
condition is (cf. ©-©) 

hm T(/x,t) = 0. (13) 

t—>--\-oo 

SOLUTION BY LAPLACE TRANSFORM 

Let F{s, t) be the Laplace transform of T{fi, t)/t with respect to fi. Then fjl2|) is equivalent 

to 

dF _ d T(0,t) _t ^. 
^~dt~di t ~ 2 ' 

dt 2s dt St ' ^ ^ 

The solution of p4j) consistent with (fT!?|) is 



00 



F{s, t) = -e^/'^ I ' e-^'/'^^ dv. (15) 

It ov sv 



F{s, t) = - T^e*'/^^ / e-'''/^^T(0, v) dv. (16) 



Equivalently, 

T(0,t) 1 
St 2s2 

Thus, in principle, T{fi,t) can be calculated from T(0,f). 

Indeed, the inverse Laplace transform can be performed at the kernel level (under the 
integral sign in ^ or (HHD) E P- 1026]: 



T(/i,t) 



fMn.V^)l{^)d. (17) 



or 



TM=T{0,t)-t -^^J,(mV^^^)T{0,v). (18) 



It ^^F^ 

A change of variable (w^ = v"^ — t^) converts (fTSj) to 

r°° mdw 



T{^,,t) = T{0,t)-t -^^J,(mw)T{0,V^i;^), (19) 
Jo \/w^ + V 

and there is a similar variant of (fTTj). 
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Example 1 



The cylinder kernel of the free massless scalar field in spatial dimension d is 

m^,x,y) = ^^i^ (20) 
where z = |x — y| is the spatial separation and c = ^{d + 1) . According to ()17|). therefore, 

T(,.(.x.y) = 2.r(cK'*^ . (21) 



From [ij, p. 425] follows 



T{m\ t, X, y) = 2i-%-Wt(t2 + z'^y'/^K,{mV^T^), (22) 

which is the correct formula^ for the cylinder kernel of the free field of mass m. All the 
solution formulas for problems with infinite fiat boundaries now follow by the method of 
images. 



Example 2 

Let d = 1 and consider an interval of length L with Dirichlet boundary conditions. For 

the massless case the solution by images can be summed in closed form p, (23) and (27)] 
with the result 

^/ N 1 sinh(7rt/L) 1 . . 

TrT O,t = , / , V • (23) 

^ ' ' 2 cosh(7rt/L) - 1 2 ^ ^ 

We apply p9|) and (jSJ and compare with the conclusions of Hays [2] about the renormalized 
total energy in the massive case (finding complete agreement). It is convenient to separate 
out the contribution of the free Green function (j^UI) (with c = 1, z = 0, and integrated over 
< X < L), because that is where all the divergences lie. That is, in both places in (fT^ 
write 



T{0,v) = — + 
irv 

a 



T(0, v)- — 
irv 



When we apply the operator — | limt_,o ^ to (fT^ we thus encounter four terms: 
• The divergent term 

(present already in empty space) is the mass-independent part of the renormalization 



of the bag constant in [2, (3.10)]. 



^ The Bender-Hays Green function |3j is the Green function of the Helmhohz equation in one higher 
dimension, [— — ^ + m?]G = 5{t)5'^{-z). In the free case this kernel is known [isi (4.25)] and is 
proportional to Kc-iim^/t"^ + z^). Differentiation and a recursion relation for the modified Bessel function 
then lead to if^ . 
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The remaining (bracket) contribution of the first term in ()19|) is the famihar massless 
Casimir energy, 

(It comes from the 0{t^) term in the Taylor expansion of ()23p . Regrettably, that 
crucial term is written in [gI, (27)] with the wrong sign.) 

The contribution of the free Green function to the integral in p9|) works out to 

rn?Lt , f mt\ , ^ ^m}Lt 
In — + (2^ 



27r V 2 / ' ' 47r 

(C = 0.577. . . , Euler's constant). The corresponding term in the regularized energy, 

m^L , (mt\ , ^ .rn?L , . 

is the mass- dependent part of the renormalization H (3.10)]. It, also, is present 
in empty space. (It includes a finite term, proportional to rn?L, which is actually 
ambiguous in the sense that the scale factor in the argument of the logarithm function 
is arbitrary.) 

The remaining (bracket) part of the integral splits into two disparate pieces. 

— The term — | in contributes^ 

_ ilini-^(l-e-™*) = -- (27) 

to the energy, in agreement with 0, (3.15)]. This constant term, associated 
with paths that reflect from the boundaries an odd number of times P, Sec. 4], 
represents the energy of interaction of the massive field with the two boundaries 
separately (i.e, it survives when L approaches infinity, and it does not contribute 
to the Casimir force). 

— What remains is the contribution of the paths that reflect an even number of 
times; it is the mass-dependent part of the true Casimir energy. In our present 
approach it equals^ 



m Ji{mw) 



4 Jo w 



smh{T{w / L) 2L 



cosh{TTW / L) — 1 nw 
u\ 2 



dw 
du 



m r°° ^ fmLu\ \ , f u\ 2] du 
In Hays's approach 0, (3.11) and (3.13)] it appears as 



2n^, s +24L ^^^^ 



This follows from (6.552.1)] and the power series of the Bessel functions. 
^ The error made by setting t = inside the integrand of this term of H19|) before differentiating is of order 
Ini, so it vanishes in the limit. 
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(since Hays's sum includes the massless Casimir energy, 1^^). The equivalence of 
()28p and ()29p is not obvious, but it can be verified by the method of H p. 427]. 
(It has also been tested numerically.) We claim no practical advantage for the 
integral, since the sum converges faster. 



PARTIAL SOLUTION BY RECURSION 

Rarely will one of the integrals (fT3jl - ()19|) be evaluatable analytically in any particular 
case. Moreover, T{0,t) often will not be available for arbitrarily large t. It is worthwhile, 
therefore, to see how much information can be obtained from the known asymptotic structure 
((Tj) if the coefficients for m = are known. 

By substituting ((Tj) into (fT^ one obtains 

U-d + s-i)^r'^--^+ t (-rf + s-i)§^t— int+ t l^t— 

s=0 s=d+l s=d+l 

s—d odd s—d odd 

OO cx) r 

= g£^t-<i+-2^ g Z^t-'i+-2^nt (30) 

s=2 ^ s=d+3 ^ 

s—d odd 

Therefore, we have the recursion relations 
for s — d odd and positive, and 

for all nonnegative integers s, the terms being set to when not defined. Generically these 
equations can be solved recursively for /^(/i) and es(/i) , respectively (the initial data /s(0) 
and 6^(0) being presumed known). Exceptions occur when —d + s — 1 = (i.e., s = (i + 1); 
then becomes a tautology and (jH^ . 

9fd+i _ Gd-i f . 

~dir~^' ^^^^ 

takes its place as the equation determining fd+i- 

Thus there is no equation to determine ed+i{fi)- The reason is that there is no way in 
this approach to impose the second boundary condition (fT!^ . so the solution must involve 
an arbitrary function. Ironically, that function turns out to be naturally identified with the 
renormalized vacuum energy (jH)), precisely the quantity of greatest physical interest. In fact, 
the only coefficients that have been completely determined by this exercise are the ones that 
are equivalent to heat-kernel coefficients 0, 0, 1^ ■ Nevertheless, the calculation clarifies the 
structure of the problem and shows that once ed+i{fi) is known (along with the mass-zero 
coefficients), all the higher cylinder- kernel coefficients are computable. 

I thank Stuart Dowker and Klaus Kirsten for comments on the first draft, and Todd 
Zapata for help with the numerical verification of "(|^H|) = ((221) "• 
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